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Abstract 

We use the fractional integrals in order to describe dynamical processes in the fractal 
media. We consider the "fractional" continuous medium model for the fractal media 
and derive the fractional generalization of the equations of balance of mass density, 
momentum density, and internal energy. The fractional generalization of Navier-Stokes 
and Euler equations are considered. We derive the equilibrium equation for fractal 
media. The sound waves in the continuous medium model for fractional media are 
considered. 
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1 Introduction 

The pore space of real media is characterized by an extremely complex and irregular geometry 
m 12 EI- Because the methods of Euclidean geometry, which ordinarily deals with regular 
sets, are purely suited for describing objects such as in nature, stochastic models are taken 
into account 0. Another possible way of describing a complex structure of the pore space 
is to use fractal theory of sets of fractional dimensionality |S]-jH]- Although, the fractal 
dimensionality does not reflect completely the geometric and dynamic properties of the 
fractal, it nevertheless permits a number of important conclusions about the behavior of 
fractal structures. For example, if it is assumed that matter with a constant density is 
distributed over the fractal, then the mass of the fractal enclosed in a volume of characteristic 
dimension R satisfles the scaling law M{R) ~ R^ , whereas for a regular n-dimensional 
Euclidean object M{R) ~ R"'. Let us assume that a network of pore channels can be treated 
on a scale i? as a stochastic fractal of dimensionality D < 3 embedded in a Euclidean space of 
dimensionality n = 3. Naturally, in real objects the fractal structure cannot be observed on 
all scales but only those for which Rp < R < Rq, where Rp is the characteristic dimensionality 
of the pore channel, and Rq is the macroscopic scale for uniformity of the investigated 
structure and processes. For example, Katz and Thompson ^01 presented experimental 
evidence indicating that the pore spaces of a set of sandstone samples are fractals in length 
extending from 10 angstrom to 100 fim. 



The natural question arises: what happens to the "laws of fluid flow" when the medium 
through which the fluid actually flows is fractional? The scientiflc reason for our interest 
in this topic is immediately obvious: the laws associated with fluid flow through fractional 
media are only beginning to be understood 

In the general case, the fractal media cannot be considered as continuous media. There 
are points and domains that are not flUed of the medium particles. These domains are the 
porous. We suggest to consider the fractal media as special (fractional) continuous media 
j24j . We use the procedure of replacement of the fractal medium with fractal mass dimension 
by some continuous medium that is described by fractional integrals. This procedure is a 
fractional generalization of Christensen approach |T^. Suggested procedure leads to the 
fractional integration and differentiation to describe fractal media. The fractional integrals 
allow us to take into account the fractality of the media. In order to describe the fractal 
medium by continuous medium model we must use the fractional integrals. 

In many problems the real fractal structure of matter can be disregarded and the medium 
can be replaced by some "fractional" continuous mathematical model. In order to describe 
the medium with non-integer mass dimension, we must use the fractional calculus. Smooth- 
ing of the microscopic characteristics over the physically inflnitesimal volume, we transform 
the initial fractal medium into "fractional" continuous model that uses the fractional inte- 
grals. The order of fractional integral is equal to the fractal mass dimension of the medium. 

More consistent approach to describe the fractal media is connected with the mathemat- 
ical deflnition the integrals on fractals. In ^2]; was proved that integrals on net of fractals 
can be approximated by fractional integrals. In ^3] , we proved that fractional integrals can 
be considered as integrals over the space with fractional dimension up to numerical factor. 
To prove we use the well-known formulas of dimensional regularizations [T^ . 

The fractional continuous models of fractal media can have a wide application. This 
is due in part to the relatively small numbers of parameters that deflne a random fractal 
medium of great complexity and rich structure. The fractional continuous model allows us to 
describe dynamics of wide class of fractal media. Fractional integrals can be used to derive 
the fractional generalization of the equations of balance for the fractal media. 

In this paper, we use the fractional integrals in order to describe dynamical processes in 
the fractal media. In section 2, we consider the "fractional" continuous medium model for 
the fractal medium. In section 3-5, we derive the fractional generalization of the equations 
of balance of mass density, momentum density, and internal energy. In section 6, the frac- 
tional generalization of Navier-Stokes and Euler equations are considered. In section 7, we 
derive the equilibrium equation for fractal media. In section 8, we consider the fractional 
generalization of Bernoulli integral. In section 9, the sound waves in the continuous medium 
model for fractional media are considered. Finally, a short conclusion is given in section 10. 

2 Fractal Media and Fractional Integrals 

The cornerstone of fractals is the meaning of dimension, speciflcally the fractal dimension. 
Fractal dimension can be best calculated by box counting method which means drawing a 
box of size R and counting the mass inside. The mass fractal dimension |22] can be easy 
measured for fractal media. 
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The properties of the fractal media hke mass obeys a power law relation 



M{R) = kR 



(D<3) 



(1) 



where M is the mass of fractal medium, is a box size (or a sphere radius), and D is a 
mass fractal dimension. Amount of mass of a medium inside a box of size R has a power 
law relation ((H). 

The power law relation (P) can be naturally derived by using the fractional integral. 
In this section, we prove that the mass fractal dimension is connected with the order of 
fractional integrals. Therefore, the fractional integrals can be used to describe fractal media 
with non-integer mass dimensions. 

Let us consider the region Wa in three-dimensional Euclidean space E^, where A is the 
midpoint of this region. The volume of the region Wa is denoted by K(WOi). If the region 
Wa is a ball with the radius Ra, then the midpoint A is a center of the ball, and the volume 
F(VFyi) = {4:/3)ttR\ . The mass of the region Wa in the fractal medium is denoted by 
Mo{Wa), where D is a mass dimension of the medium. 

The fractality of medium means that the mass of this medium in any region Wa of 
Euclidean space increase more slowly that the volume of this region. For the ball region 
of the fractal media, this property can be described by the power law (Q), where R is the 
radius of the ball Wa that is much more than the mean radius Rp of the porous sphere. 

Fractal media are called homogeneous fractal media if the power law (Q) does not depends 
on the translation and rotation of the region. The homogeneity property of the media can 
be formulated in the form: For all regions Wa and Wb of the homogeneous fractal media 
such that the volumes are equal l^(PK4) = V{Wb), we have that the mass of these regions 
are equal Md(Wa) = M£){Wb)- Note that the wide class of the fractal media satisfies the 
homogeneous property. In many cases, we can consider the porous media |15| ITB] , polymers 
[T7] . colloid agregates jUI, and aerogels [12] as homogeneous fractal media. 

To describe the fractal medium, we must use the continuous medium model such that 
the fractality and homogeneity properties can be realized in the form: 

(1) Fractality: The mass of the ball region W of fractal medium obeys a power law relation 



where D < 3 and R is the radius of the ball. In the general case, we have the scaling law 
relation 



where XW = {Xx, x G W}. 

(2) Homogeneity: The local density of homogeneous fractal medium is translation invariant 
value that have the form p(r) = po = const. 

We can realize these requirements by the fractional generalization of the equation 




(2) 



dMoiXW) = X'^dMoiW), 




(3) 
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Let us define the fractional integral in Euclidean space in the Riesz form [20] by the 
equation 

(/^p)(ro) = / pir)dVn, (4) 



w 



where dVa = C3{D,r,rQ)d^r, and 



The point Tq E W is the initial point of the fractional integral. We will use the initial points 
in the integrals are set to zero (rg = 0). The numerical factor in Eq. has this form 
in order to derive usual integral in the limit D (3 — 0). Note that the usual numerical 
factor %\D) = T{l/2)/2^n^/^T{D/2), which is used in Ref. [201, leads to 73^^(3 - 0) = 
r(l/2)/237r3/2r(3/2) in the limit D ^ (3 - 0). 

Using notations we can rewrite Eq. (jS)) in the form M^{W) = (/^p)(ro). Therefore 
the fractional generalization of this equation can be defined in the form 

Mn{W) = (/^p)(ro) = ^^^^^ jj{v)\r-v,r'dh. (5) 

If we consider the homogeneous fractal media (p(r) = po = const) and the ball region W , 
then we have 

where R = r — tq. Using the spherical coordinates, we get 

_ 7r25-^r(3/2) j _ 2^-^7rr(3/2) 

^^^^ " r(Z^/2) - DT{D/2) ' 

where R = |R|. As the result, we have M{W) ~ R^, i.e., we derive Eq. Q up to the 
numerical factor. Therefore, the fractal medium with non-integer mass dimension D can be 
described by fractional integral of order D. 

Note that the interpretation of the fractional integration is connected with fractional 
dimension jT3j. This interpretation follows from the well-known formulas for dimensional 
regularizations [Hj : 

Using Eq. (jH)), we get that the fractional integral 

f{x)dVD, 

w 

can be considered as a integral in the fractional dimension space 

^<^/^' (7) 



27r^/2r(D) 

up to the numerical factor T{D/2)/{27t°/^T{D)) 
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3 Equation of Balance of Mass Density 



The fractional integrals can be used not only to calculate the mass dimensions of fractal 
media. Fractional integration can be used to describe the dynamical processes in the fractal 
media. Using fractional integrals, we can derive the fractional generalization of dynamical 
equations ^SlEni- In this section, we derive the fractional analog of the equation of continuity 
for the fractal media. 

Let us consider the region W of the medium. The boundary of this region is denoted 
by dW. Suppose the medium in the region W has the mass dimension D. In general, the 
medium on the boundary dW has the dimension d. In the general case, the dimension d is 
not equal to 2 and is not equal to {D — 1). 

The balance of the mass density is described by the equation 

^ / p{K,t)dVD = 0, (8) 

where we use dM£){W)/dt = and the following notations: 

dVn='^-^^^p-\K\''~'dVs, dVs = d'K. 

Here, and late we use the initial points in the integrals are set to zero (fq = 0). The integral 
(jHI) is considered for the region W which moves with the medium. The field of the velocity 
is denoted by u = u(R, t). The total time derivative of the volume integral is defined by the 
equation 

^ / AdVo = [ ^dVn + [ Au^dS,. (9) 
Jw Jw '^'t Jaw 

Here, m„ is defined by Un = (u, n) = UkUk, the vector u = UkGk is a velocity field, and 
n = n^efc is a vector of normal. The surface integral for the boundary dW can be represented 
as a volume integral for the region W. To realize the representation, we derive the fractional 
generalization of the Gauss theorem (see Appendix). As the result, we have the following 
equation for the total time derivative Q of the volume integral 

jJ^AdVD = + c^\D,R)dzv{c2{d,R)Au)yVD. (10) 

For the integer dimensions d = 2 and D = 3, we have the usual equation. Let us introduce 
notations that simplify the form of equations. We use the following generalization of the 
total time derivative 

|)^ = | + c(A<i.iJK^. (11) 

where the coefficient c{D, d, R) is defined by 

oZ)-d-lp/^)/2^ 

ciD,d,R) = c^\D,R)c2{d,R) = r(3/2)r(rf/2) l^''^'"''- 
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Note that the media with integer dimensions {D = 3, d = 2) have c{D,d,R) = 1. We use 
the following generalization of the divergence 



r(3/2)r(^;/2) 

and the derivative with respect to the coordinates 

V^^-r-VD p^ dc2id,R)A _ 2^-'^-^T{D/2) d 



'12) 



Kl^-'A. (13) 



Here, DivD = ^k^k ^^"^ DivD{A\\) = V^{Auk). Note that the rule of term-by-term 
differentiation for the operator is not satisfied 

V^iAB)^AV^iB) + BV^{A). 

This operator satisfies the following rule: 

V^iAB) = AV^iB) + ciD, d, R)B^^. (14) 

Note that Vf (1) 7^ 0, and we have 

V^(l) = c(D, d, R){d - 2)xk/R\ 

Using these notations, we rewrite Eq. ()1()|1 for the total time derivative of the integral in an 
equivalent form 

S i^"^- = + D,M<^)yv,. (15) 

To derive the fractional generalization of the equation of continuity [TB], we consider 
A = p(R, t) in Eq. (fT3j) and the equation of balance of mass (jH)). Substituting A = p(R, t) 
in Eq. (fT3j) . we get 

Therefore, equation of balance of mass (jH)) has the form 



This equation is satisfied for all regions W. Therefore, we have the fractional equation of 
continuity 

[^)^p + pDtVD{u) = 0. (17) 
The fractional generalization of equation of continuity has the form 



Using Eqs. (jHTIjl and (ffn|) of Appendix, the equation of continuity can rewritten in an 
equivalent form 

Using 

we have the equation of continuity in the form 

^ + c{D, d, R)u^ + ciD, d, R)p(dtv{u) + id- 2)^1^) = 0, (19) 

where we use the following notation (u, R) = UkXk- 

For the homogeneous media, we have p = const and the equation of continuity leads us 
to the equation 

, , , (rf-2)(R,u) 
diviu) + ^ ^^1^ ' ^ = 0. 

Therefore, we get the non-solenoidal field of the velocity {divu = du/dH. ^ 0). 

In addition to mass density p, the continuity equation includes the density of momentum 
pu. To obtain the equation for the density of momentum, we consider the mass force and 
surface force. 



4 Equation of Balance of Momentum Density 

Let the force f = fk^k be a function of the space-time point (R, t). The force F^, that acts 
on the mass M£,{W) of the medium region W, is defined by 

F*^= [ p(R,t)fiR,t)dVD. (20) 
Jw 

The force F*^, that acts on the surface of the boundary dW of continuous medium region 
W, is defined by 

F^= / p„(R,t)dS'rf, (21) 

JdW 

where p = p(R, t) is a density of the surface force, and p„ = Pki^k^k- Here, n = Uk^k is the 
vector of normal. 

Let P be a momentum of the medium mass that is situated in the region W . If the 
mass dMrtiW) = p(R,t)dVD moves with the velocity u, than the momentum of this mass 
is dP = dMD{W)u = pudVo- The momentum P of the mass of the region W is defined by 
the equation 

piR,t)u{R,t)dVD. (22) 



'w 

The equation of balance of density of momentum 



^ = F^^ + F^. (23) 



Substituting Eqs. - into Eq. (f^ . we get the balance equation in the form 

/ p(R,t)u{R,t)dVD= [ p(R,t)f{R,t)dVD+ [ Pn{'R,t)dSd. (24) 

Jw Jw JdW 

Using Eq. (ffT|) of Appendix, the surface integral can be represented as the following volume 
integral 

/ PndS,= [ C2{d, R)p^dS2 = [ ^^^^^^^Mc-\D,R)dVD= [ Vf'pdVD. 
Jaw JdW Jw (^^i Jw 

The balance equation of momentum density can be written for the components of vectors 
u = MfcCfc, f = fkSk, and pi = PkiSk, in the form 



/ pukdVo 
w 



dt 



[ (pfk + Vf'pki)dVD. 
Jw ^ ' 



The relation for the derivative of the volume integral with respect to time has the form ()15|1 . 
Using this relation for A = puk, we get the equation for the total time derivative of the 
integral 

— J pukdVD = j {Jy—^J^pUk) + {pUk)DivD{\y)^dVD. (25) 
Therefore, we can rewrite the balance momentum density in the form 
f d 

I iidi) D^'^^''^ ^'^''^ Dwd{u) - pfk - Vfpki^dVo = 0. 
This equation is satisfied for all regions W. Therefore, 

(^)^(pMfc) + (pUk) DivDiu) - pfk - Vfpki = 0. 
Using the rule of the term-by-term differentiation with respect to time 

we get the following form of the equations 

+ Uk(^(^^^ ^p + p Divoin)^ - pfk - ^fpki = 0. 



/ d\ / / d 

P\dt)n''' 

Using the continuity equation, we reduce the fractional equation of balance of density of 
momentum to the form 

p(^)^Wfc = p/fc + Vfpfc;. (26) 
These equations can be called the equation of balance of momentum of fractal medium. 
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5 Equation of Balance of Energy Density 

In the general case, the internal energy for the inhomogeneous medium is a function of the 
space-time point (R, t): e = e(R, t). The internal energy dE of the mass dM^^W) is equal 
to dE = e{R,t)p{Ii,t)dVD- The internal energy of the mass of the region W is defined by 
the equation 

E^ p{R,t)e{R,t)dVD. 
Jw 

The kinetic energy dT of the mass dMoiW) — pdVo, which moves with the velocity u = 
u(R, t), is equal to 

2 2 

dT = dMoY = ^d^D. 
The kinetic energy of the mass of the region W is 

T^ [ P^dVo. 
Jw 2 



The total energy is a sum of the kinetic and internal energies 

U = T + E = j^p{^ + eyVD. 
The change of the total energy is defined by 

U{t2) - U{ti) ^Am + As + Qs: 

where Am is the work of mass forces; Am is the work of surface forces; Qs is the heat that 
are influx into the region. 

The mass dMo{W) = pdVo is subjected to force ipdVo- The work of this force is 
{vL^i)pdVDdt, where (u, f) = Ukfk- The work of the mass forces for the region W and time 
interval [ti;t2] is defined by the following equation 



t2 r 

dt I {n,t)pdVD. 
ti Jw 



The surface element dS^ is subjected to force PudS^- The work of this force is (p„, u)dSddt. 
The work of the surface forces for the region W and time interval [ti;t2\ is defined by the 
following equation 

As^ dt (u, Pn)dSd. 

Jti Jaw 

The heat that are influx into the region W through the surface dW is defined by 

Qs^ [ dt f qndSd, 
Jti Jaw 

where Qn — (n, q) — UkQk is the density of heat flow. Here, n is the vector of normal. 
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The velocity of the total energy change is equal to the sum of power of mass force and 
the power of surface forces, and the energy flow from through the surface: 



dt 



p[\ + <^)dVD= I (u, ^)pdVD + / (u, ip^)dSd + / qndSd. (27) 
^ ^ Jw Jaw Jaw 



'W 

Using Eq. (fT3j) for A = p(u^/2 + e), we can rewrite left-hand side of Eq. ^I7\i in the form 



Using the equation of continuity, we get 



,2 



1//(t-)--I(K1)„(t-))* 



The surface integrals in the right-hand side of Eq. ()27|) can be represented as volume integrals 

/ {n,pn)dSd= I V^{vun)dVD, (29) 
Jaw Jw 



and 



qndSd= / V^qkdVo. (30) 
aw Jw 



Substituting Eqs. (j2Hl) - (EOI) in Eq. (jZ3), we get 

7 7 /» 

pu (-) + p (-) ^e(R, t)) dVn = ((u, f )p + VP(p,, u) + Vf g,) rfVz,. (31) 



'w 

These equations can be rewritten in an equivalent form 

f d d C 

j {'^''^{jt) D^^^ ^{jt) d'^)'^^'^ ^ j {p^>^f>^ + ^?^PkiUk) + '^kqk)dVD. (32) 

Let us use equations of balance of momentum (pUj). Multiplying both sides of these equations 
on the components Uk of vector u and summing with respect to k from 1 to 3, we get the 
equation 

pUk (y^^ = pukfk + UkVfpkh (33) 
Substituting Eq. (jHSI) in Eq. we get 

j (pUkfk + UkW^Pki + p(^j^ J)dVD = j [pUkfk + ^?{PkiUk) + ^kqk)dVD. (34) 
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Using Eq. (fT^ in the form 

OXi 

we obtain 



7ii^-V?»)dVD = 0. (35) 

This equation is satisfied for all regions W . Therefore, we get the fractional equation of 
balance of density of energy in the form 



p(^)^e = c{D,d,R)pui^^+V^qu. (36) 

6 Fractional Navier-Stokes and Euler Equations 

In Sections 3-5, we derive the fractional generalizations of the balance equations for fractal 
media. As the result we have the fractional hydrodynamic equations in the form: 

(1) The fractional equation of continuity 

(2) The fractional equation of balance of density of momentum 

p{j^Uk = pfk + ^?Pki. (38) 

(3) The fractional equation of balance of density of energy 

e = c{D,d,R)pki^ + V^qk. (39) 



\dt/ D OXi 

Here, we mean the sum on the repeated index k and / from 1 to 3. We use the following 
notations 

VPA = a(D,d)R^-''^(R''~^A). (40) 
oxk V / 

iii)u-m^ - Ft + 

where a{D, d), and c{D, d, R) are defined by the equations 



c{D,d,R) = a{D,d)R'^+^-^, a{D,d) = R 



r(3/2)r(rf/2) ' 



^2 



The equations of balance of density of mass, density of momentum and density of internal 
energy makes up a set of five equations, which are not closed. These equations, in addition 
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to the hydrodynamic fields p(R, t), u(R,t), e(R, t), include also the tensor of viscous stress 
Pki(R,t) and the vector of thermal fiux qk(R,t). 

Let us start with the definition of tensor pki = pki(R,t). According to Newton's law, 
the force of viscous friction is proportional to the relative velocity of motion of medium 
layers (that is to the gradient of the relevant component of velocity). We further assume 
that tensor Pki(R., t) is symmetrical, and characterizes the dissipation due to viscous friction. 
The most general form of tensor of viscous stress, which satisfies the above requirements, is 
determined by two constants (/i and ^) which can be chosen so that 

fduk dui 2 dum\ , dum 
Pki = -pOki + /i ^ h ^ — + iOki^ — . (42) 

\OXi OXk O OXm-' OXm 

This expression contains a second coefficient of viscosity ^, called the coefficient of inter- 
nal viscosity because it refiects the existence of internal structure of particles. In case of 
structureless particles ^ = 0. 

The definition of the vector of heat fiux = Qk^R-jt), is based on the empirical Fourier 

law 

dT 

qk = -A—, (43) 

OXk 

where T = T(R, t) is the field of temperature. The value of heat conductivity A can be 
found experimentally. 

Now we have a closed set of Eqs. (gH), and (gSl) for fields p(R,t), Uk(R,t), 

T(R, t) - a set of fractional equations of hydrodynamics. The fractional generalization of the 
equations of theory of elasticity for solids can be obtained in a similar way. 

Let us consider the special cases of the set of Eqs. (jHTfl - flH^ . (j^ . and PHj) . 

1. Let us consider the fiuids that are defined by 

Pki = -phu qk = 0, 

where p = p(R, t) is the pressure. The fractional hydrodynamic equations for these fiuids 
have the following form 

'4) P = -P^kUk. (44) 



dtJ D 
d\ „ 1 



Uk = fk--V'^p. (45) 
at/ D p 

These equations are the fractional generalization of the Euler equations. 

2. Let the coefficients p, ^ and A are constants. If we consider homogeneous viscous 
fluid, then we have 

For the fractal media, we have the non-solenoidal field of the velocity {div{u) ^ 0), that 
satisfies the relation 

. ^ dum {2-d)xkUk 
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Using Qk = and Eq. , we get the fractional generalization of Navier-Stokes equations 
in the form 

+ c{D, d, R)u,—) = P/^ - V. P + /. V, — + /. V, — + - -/.) V, — . (48) 

Equations (jlTj) and (jlHj) form the system of 4 equations for 4 fields Mi(R, t), M2(R, t), M3(R, t), 
and p(R, t). Note that Eq. (PHj) can be rewritten in an equivalent form 

fduk /n ^ DA ^^'^A 



= pf,-c{D,d,R)-^^+fic{D,d,R)-^^ + {^ + ^)c{D,d,R)-^^ (49) 
Here, we use the following notations: 

Lk{D, d, R,u)=f,[^ + - ^^^(1) + - r^^,"^^^^^^ 

where Vf(l) = c^^{D,R)dc2{d,R)/dxk. lic{D,d,R) = 1 and Lk{D,d,R,u) = 0, then Eq. 
has the usual form of the Navier-Stokes equations. 



7 Fractional Equilibrium Equation 

The equilibrium state of media means that we have the conditions 

_ '^^ - 

dt ' dxk 

for the hydrodynamic fields A = {p,Uk,e}. In this case, the fractional hydrodynamic equa- 
tions have the form 

pfk + ^?Pki = 0, V^gfc = 0. 

Using dui/dxk = 0, we get that the tensor pki has the form pki = —pSki- Using the Fourier 
law, we have the following system of equations: 

fk = ^V^p, (50) 
8T 

Vf(A^) = 0. (51) 

which are the fractional generalization of the equilibrium equations for the media. Eq. (jKn|l 
can be rewritten in an equivalent form 

d{c2id,R)p) fr> mf 
-^^ = pcsiD,R)fk. 
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Let us consider the homogeneous medium with the density p{x) = const. In this case, we 
have the equation 

^ / n d{c2id, R)p/ po) 

C3{-l-^,-n-)jk = ^ • 

If the force fk is a potential force such that c^i^D, R)fk = —dU/ dxk, then we get the fractional 
generalization of equilibrium equation in the form 

C2{d, R)p + pqU = const. (52) 



8 Fractional Bernoulli Integral 



Let us consider the equation of balance of momentum density with the tensor pki = —pSki. 
Using the relation 



dt)D-2='''\dt)n'''^ 
and Eq. (jHS|) . we get 

d 1 

f:^) -77 = ^kfk- -UkVkP- (53) 
\dt/ D 2 p 

If the potential energy U and pressure p is time-independent fields {dU/dt = 0, dp/dt = 0), 
then we can use the following relation: 

Let us consider the non-potential force that is described by the equation 

fk = -c{D, d, R)dU/dxk. (55) 

If D = 3 and d = 2, then this force is potential. Using Eqs. ()54|1 and ()55j) . we can rewrite 
Eq. (j^Hj) in the form 

7 2 

^('^ + U{D,d) + P{d)) =0, 



dt\2 

where the function P is defined by the usual relation 

dic2id, R)p) 



P{d)- 

As the result we have that the integral 



Po 



C2{d, R)p 



3 2 



Y^^ + U + P{d) = const. 



k=l 

can be considered as a fractional generalization of Bernoulli integral for fractal media. If the 
forces fk are potential, then the fractional analog of the Bernoulli integral does not exists. 
If the density is described by 

p = poC2\d,R) = p,^-^!^R^-\ (56) 
14 



then we have the fractional generahzation of Bernoulh integral in the following form 

h poU{D, d) + C2(rf, R)p = const. 

If Uk = 0, then we get equilibrium Eq. (jK^ for nonpotential force and the density (fHHj) . 

9 Sound Waves in Fractal Media 

Let us consider the motion of the medium with small perturbations. Fractional equations of 
motion (jl^ and have the form 

^ + c{D, d, R)ui^^ = -pV^Uk, (57) 
— + c{D,d,R}ui— = fk~ -VkP. 58 

Ot OXi p 

Let us consider the small perturbation 

p = po + p', P = Po+ p', Uk = Mfc, (59) 

where p' « po, and p' « po. Here, po and po describe the steady state that is defined by 
the conditions 

^ = ^ = — = — = 

dt ' dxk ' dt ' dxk 

Suppose that fk = 0. Substituting (jKHjl in Eqs. (|F7jl and (jSHI), we derive the following 
equations for the first order of the perturbation: 

% = -PoV^'„ (60) 

^ = --Vfp'. (61) 
dt Po 

These equations are equations of motion for the small perturbations. To derive the indepen- 
dent equations for perturbations, we consider the partial derivative of Eq. (lUU]) with respect 
to time: 

r)ii' 

' p^r-^- (62) 



dt^ ^ ' dt 
Substituting Eq. in this equation, we get 

rf n' 



dt^ 



V-V-p'. (63) 



If we consider the adiabatic processes, we can use equation p = p{p, s). For the first order 
of perturbation we have the equation 



/ 2 / 1 1 dp 



15 



As the result, we obtain the following fractional generalization of the wave equations: 

% - ^ Vf p' = 0, (64) 



v'^'^^\:p' = 0. (65) 

Let us consider the simple example of the fractional wave equations. If we consider one 
dimensional case (n = 1), where D < 1 and C2 = 1, then we get the following equation for 
pressure 

^'("•^)^-"K^'(«-^>f)=°- («) 

where the coefficient Ci{D,x) is defined by the relation 



l-D-l 

Ci{D,x) 



\x\ 



V{D)- 

10 Solution of the Fractional Wave Equation 

The fractional generalization of wave equation has the following form 

-<---'l^-|^("^-(--)|)' 

where Ci(D, x) > 0. Let us consider the region < x < / and the following conditions: 

u{x,0) = f{x), {dtu){x,0) = g{x), 

u{0,t) = 0, u{l,t)=0. 
The solution of Eq. ()67|) has the form 

oo 

u{x,t) = ^(^fnCos{Xnt) + ^^sm( A„t) j ?/„ (x) . 



n=l 



Here, and gn are the Fourier coefficients for the functions f{x) and g{x) that are defined 
by the equations 

fn = \ \yn\\~^ / f{x)yn{x)dlD = \ \yn\\'^ / Ci{D , x) f {x)yn{x)dx, 

Jo Jo 

gn=\\yn\\''^ f{x)yn{x)dlD=\\yn\t'^ j Ci{D , x) f {x)yn{x)dx, 

Jo Jo 

\\yn\\^ = / yl{x)dlD= / ci{D,x)y'^{x)dx, 
Jo Jo 
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where dlo — ci{D,x)dli, and dli — dx. Note that the eigenfunctions yn{x) satisfy the 
following condition 

rl 

yn{x)ym{x)dlD = Snm- 







The eigenvalues A„ and the eigenfunctions yn{x) are defined as solutions of the equation 

v^[ci{D,x)y',l + X''ci{D,x)y^O, ?/(0) = 0, y{l) ^ 0. 

This equation can be rewritten in an equivalent form 

v^xy'Ux) + (D- l)y'^{x) + \^xy{x) = 0. 

The solution of this equation has the form 

y{x) = C^x'-^'/^MXx/v) + C2X^-''/^Y,{Xx/v), 

where u — |1 — D/2\. Here, Ju{x) are the Bessel functions of the first kind, and Yi,{x) are 
the Bessel functions of the second kind. 

As an example, we consider the case that is defined by 

Z = v^l, 0<x<l, f{x)^x{l-x), g{x)^0. 

The usual wave has D = 1 and the solution 

4(1 — {—l)"')sin{Trnx)cos{Trnt) 



oo 

U{X, '] 

n=l 



The approximate solution for the usual wave with D — 1 that has the form 

, , 4(1 — (—iy'')sin('Knx)cos(7rnt) 



n=l 



is shown in Fig. 1 for < t < 3 and velocity v = 1. 

li D — 1/2, then we have the fractal medium wave with 

Vnix) = :^^x'/^Js/^{^/2X^x/2). 

The eigenvalues A„ are the zeros of the Bessel function 

An : J3/4(\/2A„/2) = 0. 

For example, 

Ai - 4.937, As ~ 9.482, A3 ~ 13.862, A4 ~ 18.310, A5 ~ 22.756. 
The approximate values of the eigenfunctions 



fn — 



y j\''/\l-x)J^,^{V2Xj2)dx, 



V{D) 
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are following 



/i ~ 1.376, /2 



0.451, /3 ~ 0.416, U 



0.248, /s ~ 0.243. 



The solution of the fractional equation 



oo 



U 



{X,t) = fnCOs{\nt)J3/iiV2\nX/2). 



n=l 



The approximate solution for the fractal media wave with D = 1/2 that has the form 



is shown in Fig. 2 for the velocity v = 1. 

11 Conclusion 

The fractional continuous models of fractal media can have a wide application. This is due 
in part to the relatively small numbers of parameters that define a random fractal medium of 
great complexity and rich structure. The fractional continuous model allows us to describe 
dynamics for wide class fractal media. 

In many problems the real fractal structure of matter can be disregarded and the medium 
can be replaced by some "fractional" continuous mathematical model. To describe the 
medium with non-integer mass dimension, we must use the fractional calculus. Smoothing 
of the microscopic characteristics over the physically infinitesimal volume transform the 
initial fractal medium into "fractional" continuous model that uses the fractional integrals. 
The order of fractional integral is equal to the fractal mass dimension of the medium. 

The experimental research of the hydrodynamics of fractal media can be realized by 
introducing a neutral indicator, its distribution. It allows one to obtain much information 
on the motion and mixing of fluid. The fractal nature of damage and porosity has been 
experimentally detected over a wide range of scales. 

Note that the fractional hydrodynamic equations for fractal media can be derived from 
the fractional generalization of the Bogoliubov equations that are suggested in • 

I would like to thank Prof. G.M. Zaslavsky for very useful discussions. 



Appendix: Fractional Gauss Theorem 

To realize the representation, we derive the fractional generalization of the Gauss theorem 



10 




n=l 
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where m„ is defined by m„ = (u, n) = UkUk, the vector u = Uk^k is a velocity field, and 
n = UkCk is a vector of normal. Here, 

, , , d{Au) d{Auk) 

Here, and later we mean the sum on the repeated index k and / from 1 to 3. 
Using the relation 

dSd = C2{d, R)dS2, C2{d, R) = :^L_^\Rf-\ (69) 

we get 

/ AundSd= / C2{d, R)AundS2. 
Jaw Jaw 

Note that we have C2(2, R) = 1 for the d = 2. Using the usual Gauss theorem we get 

C2{d, R)AundS2 = / div{c2{d, R)Au)dV3. 
Jw 

The relation 

dVo = C3{D, R)dV^, c^iD, R) = 

in the form dVs = c^^(D, R)dVD allows us to derive the fractional generalization of the Gauss 
theorem: 

(71) 



aw 



AundSd= / c^^{D,R)div{c2{d,R)Au)dVD- 
Jw 
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